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ABSTRACT 

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) 
IRAS 17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited (v 2 =l) 7=3-2 and 4-3 HCN. The 
emission is emerging from buried, compact (r <17-70 pc) nuclei that have very high implied mid-infrared surface brightness >5 X 10^^ L© kpc"^. 
These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, vi^l, lines 
of HCN are excited by intense 14 pm mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG 
nuclei when H 2 column densities exceed 10^"^ cm“^. It is clear that these lines open up a new interesting avenue to gain access to the most obscured 
AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows 
for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state (y=0), 7=3-2 and 4-3 
rotational lines of HCN and HCO^ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. 
The HCN and HCO^ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form 
of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we 
discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear 
growth, before the phase of strong feedback. 

Key words, galaxies: evolution — galaxies: individual: Arp220, IRAS 17208-0014, IC860, Zw049.057 — galaxies: active — galaxies: nuclei — 
galaxies: ISM — ISM: molecules 


1. Introduction 


Luminous (Lir= 10^^ - 10^^ L©) and ultraluminous (Ljr^IO^^ 
Lq) infrared galaxies ((U)LIRGS) are powered by either bursts 
of star formation or AGNs (active galactic nuclei - accreting 
supermassive black holes (SMBHs)) and a re fundamental to 
galaxy mass assembly over cosm ic time (e.g. lElbaz & Cesarskvl 
120031: ISanders & Mirabellfl^ . Some (U)LIRGS have deeply 
embedded nuclei that harbour a very active evolutionary stage 
of AGNs and starbursts, often with signatur e s of outflow¬ 
ing and infalling gas (e.g. iBanerii et al.l l20lf 


_^ _ lEabiai] 119991: 

Gonzalez-Alfonso et al.l 120121: ISilk & Reesl 119981: ISpoon et al.l 
20131) . These nuclei are key to understanding nuclear growth and 


feedback mechanisms and studying them is essential for a com¬ 
plete AGN and starburst census, for constraining orientation- 
based unification mod els and fully understanding how galaxies 
grow and evolve (e.g. iBrightman & Uedal 120121: iMerloni et al.l 
[ 2011 . 


However, probing inside the optically thick layers of 
the most enshrouded galaxy nuclei is an observational chal¬ 
lenge. X-ray and mid-infrared (mid-IR) surveys are designed 


Based on observations carried out with the IRAM Plateau de 
Bure and ALMA Interferometers. IRAM is supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). ALMA is a partnership of 
ESO (representing its member states), NSF (USA), and NINS (Japan), 
together with NRC (Canada) and NSC and ASIAA (Taiwan), in coop¬ 
eration with the Republic of Chile. The Joint ALMA Observatory is 
operated by ESO, AUI/NRAO, and NAOJ. 


to find embedded nuclear activity, but when Ay 


1000 


magnitudes (/V(H 2 )> 10 cm“^), the emission becomes ex¬ 

tremely attenuated even at these wavelengths (e.g.iTTeister et af] 
I 2 OIOI: ILusso et ^120131: ISpoon et aP 120021: iRoche et al.ll2015l) . 
Continuum ima ging at the more tra nsparent radio- (e.g . 
iParra et aT]|2010l) and mm/sub-mm (e.g. ISakamoto et all 120081) 
wavelengths can be used to probe luminosity density and to iden¬ 
tify power sources, although interpretations are sometimes com¬ 
plicated by surrounding star formation, free-free absorption, and 
other optical depth effects. Rotational lines (in the vibrational 
ground state y=0) of polar molecules, such as HCN and HCO^, 
are common tracers of dense (n > 10"^ cm“^) gas and are of¬ 
ten used t o probe the nuclear mass, dynamics, and astroche m- 
istry (e.g. iGao & Solom(^l2004 iMeijerink & Spaansll2005l) of 
obscured galaxy nuclei. However, large optical depths and the 
low-level energies of the transitions (see Eig.[T]) make them less 
well suited to studying deeply embedded, hot, and compact star- 
bursts or obscured AGNs. 


In contrast, rotational transitions of vibrationally excited 
(y 2 =l) HCN occur b etween energy levels exceeding 1000 K 
above ground (EigdJ (IZiurvs & Turnej[r986h . These transitions 
are usually radiatively excited by the intense mid-IR emission 
provided by the presence of hot dust and high column densi¬ 
ties (see Sect. [2| and Appendix |B]). Therefore, the dust and gas 
obscuration impeding the use of other tracers is a prerequisite 
for the excitation of vibrational HCN. Since the rotational tran¬ 
sitions of HCN occur in the mm and sub-mm regimes the line 
emission can better penetrate the optically thick layers of dust. 
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Fig. 1. HCN energy diagram with the HCN y=0 and y 2 =l 
J=4,3,2 levels marked. The split of the y 2 =l rotational lev¬ 
els into e and / is also shown. The energy levels are not to 
scale which is illustrated by the dashed y-axis. The transitions 
are indicated with arrows. The values for the energy levels are 
from the Cologne Database for Molecule Spectroscopy (CDMS) 
(http s ://w w w. astro. uni-koeln. de/cdms) 


so it functions as an unattenuated tracer of deeply buried high 
mid-infrared surface brightness nuclei. 


Extragalactic vibrationally excited HCN h a s bee n detected 
in the LIRG NGC 4418 dSakamoto et all 1201(1) and the 
ULIRGs Arp220 (ISalter et all l2008h . IRAS 20551-4250 
dlmanishi & Nakanishil l2013l) . and Mrk231 (lAalto et alJ 
l201.5h . In NGC4418 and Arp220, vibrational lines of HC 3 N 
have also been used to stud y the p roperties of the buried 


nucleus (ICostagliola & Aalt^ 120101: Sakamoto et al.l 120101: 


Martin et all 120111: ICostagliola et al.l 20131: ISakamoto et^ 


_ igl; 

20131: ICostagliola et al.|[2015l) . Vibrational temperatures ranging 


from 200 to 400 K have been estimated for HCN and HC 3 N 
in NGC4418, Arp220 and Mrk231. Vibrationally excited HCN 
has been detected in the circum-nuc lear disk at the centre of 
the Milky Way bv iMills et al.l (l2013l) . but line luminosities are 
significantly lower (relative to those of y=0 HCN) compared to 
those in the galaxies studied here. 


In this paper we present IRAM Plateau de Bure A-array ob¬ 
servations of / = 3 ^ 2, y-0 HCN, HCO+ and y 2 =l HCN 
emission in the LIRGs IC860 and Zw049.057. We also present 
Atacama Large Millimeter Array (ALMA) observations of the 
corresponding / = 4 ^ 3 lines in the ULIRGs Arp220 and 
IRAS 17208-0014. These new observations double the num¬ 
ber of galaxies with detected vibrationally excited HCN emis¬ 
sion and reveal in unprecedented detail emission emerging from 
highly compact structures behind a temperature gradient. We can 
show, for the first time, that the ground vibrational state HCN 
and HCO^ line emission becomes self- and continuum absorbed 
and fails to trace the inner, nuclear regions of these galaxies. 


2. Vibrationally excited HCN 

The HCN molecule has degenerate bending modes in the IR and 
may absorb IR-photons to its first vibrational state (Lig.[T]). The 


first excited bending state of HCN is doubly degenerate, so that 
each rotational level y 2 =l / is split into a closely spaced pair of 
levels (/ and e). The J=3-2 and 4-3 y 2 =l lines are rotational 
transitions within the vibrationally excited state that occur near 
the frequency of the corresponding y=0 lines. Lor example, the 
y 2 =l^ 3-2 line is shifted by only 20 MHz from the y=0 HCN 3- 
2 line and therefore appears close to the centre of the HCN line 
for extragalactic sources while the y 2 =l/ line has rest frequency 
v=267.1993 GHz, shifted by 358 MHz (-1-400 kms“^) from the 
HCO^ 3-2 line. The y 2 =l^ and If line emission always have 
similar intensity (except for cases of population inversion). 

We will hereafter refer to the HCN y=0 lines as HCN and the 
vibrationally excited HCN y 2 = If lines as HCN-VIB. 

2.1. Mid-IR excitation of the HCN-VIB lines 

The energy above ground of the y 2 =l state starts at r£/j^=1024 
K and the critical density of the HCN-VIB lines is large 
enough (> 10 ^^ cm~^) that collisional e xcitation is unlikely 
(IZiurvs & Turner! Il986l: iMills et aP |2013|) . Instead, the vibra¬ 
tional ladder can be excited through the absorption of mid- 
IR 14yum continuum emission. There is no sharp threshold, 
but at mid-IR background temperatures of TB(l4/jm)^100 K 
the pumping is efficient enough to start t o excite the line 
dCarroll & Goldsmithll1981l:lAalto et alJl2nn7h . 

Lor a thorough investigation of the vibrationally excited 
HCN a multi-line, non-LTE analysis is required. Lor exam¬ 
ple, if the buried source is very hot, compact and luminous it 
may be necessary to also consider higher vibrational modes (yi 
at 3.02jum and y^ at 4 .77jum) for the excitation of HCN-VIB 
(ICernicharo et akOOl ll) . In a forthcoming paper we will present 
non-LTE models of HCN-VIB emission excited by buried IR 
sources (Gonzalez-Alfonso et al. (in prep.)). 

However, for the purpose of this discussion we will assume 
that HCN-VIB is excited primarly by 14/im continuum with a 
brightness temperature Tb ~ 100 K. 


3. Observations 

3.1. Sample galaxies 

The studied galaxies (see Tab. [T]) are part of a sample where 
we find evidence of compact nuclear activity wi th faint 158 um 
[C II1 line emission relative to LIR continuum (iMalhotra et al.l 
1 19971 iDiaz-Santos et al.l l2013l) and/or th at show prominent 
14yum HCN absorption (iLahuis et al.ll2007l) . The [C II] deficit 
correlates with de gree of compactness of t he IR emitting source 
and the L^irIMh^ (iDiaz-Santos et al.|[2Q13l) . 

The studied galaxies also have unusually bright HC 3 N line 
emission which we prev iously have suggested is due to a com¬ 
pact buried hot nucleus (I Aalto et al .1120071: ICostagliola & Aalt^ 
I 2 OIQI: ICostagliola et all l201ll) . Note that we only include the 
western nucleus of the ULIRG merger Arp220 in this study 
(hereafter Arp220W) which is the brighter of the two nuclei. The 
eastern nucleus is discussed in Martin et al. (in prep.) 

3.2. IRAM Plateau de Bure 

The 1 mm observations were carried out with the 6 -element ar¬ 
ray in March 2014 in A-array. The receivers were tuned to a fre¬ 
quency of 263.0 GHz, to include both the redshifted HCN and 
HCO+ J = 3-2 lines (y=265.886 GHz and 267.558 GHz rest 
frequency, respectively). We used the WideX correlator provid- 
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Table 1. Sample galaxies^ 


Name 

R.A. 

(J2000) 

Dec. 

(J2000) 

Log Lfir 
[Lo] 

D 

[Mpc] 

Type 

Telescope 

Beam 

(") 


Sensitivity* 
(mJy / channel) 

IC860 

13:15:03.05 

+24:37:08.0 

11.17 

59 

Sab 

IRAM PdBI 

0."50 X 0.'' 

’23 

1.2 mJy 

Zw049.057 

15:13:13.10 

+07:13:32.0 

11.27 

59 

Sab 

IRAM PdBI 

0."55 X 0." 

’24 

1.6 mJy 

Arp220W 

15:34:57.27 

+23:30:10.5 

12.03 

80 

merger 

ALMA 

0/'44 X 0.'^ 

’38 

1.9 mJy 

IRAS 17208-0014 

17:23:21.94 

-00:17:01.0 

12.39 

176 

merger 

ALMA 

0."49 X 0." 

’34 

1.6 mJy 


^Lir and D from I Sanders et alJ j2003h (where D is the angular distance {Da) ) but for Arp220W we attempt to estimate the lumi¬ 
nosity emerging from the western nucleus only from the continuum study of I Sakamoto et alJ (l2008h . For channel widths of 33 
km s“\lC860 and Zw049.057) and 20 km s"\Arp220W and IRAS 17208-0014). 


ing a broad frequency range of 3.6 GHz and 2 MHz spectral res¬ 
olution. After calibration within the GILDAS reduction package, 
the visibility set was converted into FITS format, and imported 
into the AIPS package for further imaging. We adopted natural 
weighting during the CLEAN deconvolution. Table [T] gives the 
beam sizes and sensitivities of the final data. The details of the 
observations are: 

IC860: The bandpass of the individual antennas was derived 
from the bright quasar 3C84 and the absolute fiux calibration 
was set on MWC349. The quasars JI310-f323 (~ 0.7 Jy at I mm) 
and JI328-1-307 (~ 0.2 Jy) were observed regularly for complex 
gain calibration every 25 minutes. 

Zw049.057: The bandpass and fiux calibrations were derived 
from 3C279 and MWC349, respectively. The gain calibration 
was derived from observations of JI502-I-I06 (~ 0.4 Jy) and 
JI456-I-044 (~ 0.2 Jy) with a similar scheme as for IC860. 

3.3. ALMA 

The final sensitivity and beam sizes of the observations are given 
in Table [Hand the details are discussed below. 

IRAS17208-0014: The 0.8 mm observations were carried out 
with ALMA (with 34 antennas in the array) in May 2014. 
The correlator was set up to cover two contiguous bands of 
1.875 GHz in each of the lower (centered at ~341.9 GHz) and 
upper (centered at -352.1 GHz) sideband, with 960 spectral 
channels (1.95 MHz wide each) per band. The upper side band 
includes the HCN / = 4-3v = 0 and v = 1/ lines. The bandpass 
of the individual antennas was derived from the bright quasar 
JI733-I304. The quasar JI75I-I-096 (~ 1.7 Jy) was observed 
regularly for complex gain calibration. Its fiux was extracted 
from the ALMA fiux-calibrator database for absolute fiux cal¬ 
ibration. After calibration within the CASA reduction package, 
the visibility set was imported into the AIPS package for further 
imaging. 


Arp220: The 0.8 mm observations were carried out with 
ALMA in April and May 2014 (with 34 and 30 usable antennas, 
respectively). Spectral features presented here were observed in 
the USB with 2 x 1.875 GHz spectral windows centered in 347.8 
and 349.6 GHz. The bandpass of the individual antennas was de¬ 
rived from the quasar JI550-f0527. The fiux calibration was set 


on Titan, where intrinsic atmospheric lines were masked. The 
quasar JI5I6-f 1932 (~ 0.4 Jy at 0.9 mm) were observed regu¬ 
larly for complex gain calibration every 5 minutes. Calibration, 
deconvolution and imaging were carried out in the CASA reduc¬ 
tion package. 

4. Results 

Line fiuxes, luminosities and source size are presented in Tab.|2l 
Nuclear (central beam) spectra in Fig.|2|and Fig.[3l 

4.1. HCNandHCO^ 

HCN and HCO^ source diameters (FWHM) range between 
100 and 400 pc (Tab. O and the line profiles have a striking 
double-peaked shape with a minimum near systemic velocity 
(see Figs. (Hand O. 

4.1.1. Line shapes 

In Sect. 15.11 we present strong evidence that the double-peaked 
HCN and HCO^ spectra are caused by deep self- and contin¬ 
uum absorption at systemic velocity. Therefore, FWHM widths 
cannot be fitted for these lines. In IC860, the shapes of the HCN 
and HCO^ spectra appear as reversed P-Cygni profiles - strongly 
asymmetric with the red part of the line profile showing stronger 
absorption. This is particularly striking for HCO^ where the con¬ 
tinuum is also absorbed on the red side. The deepest absorption 
feature occurs at a projected velocity of 70 kms“^ The emis¬ 
sion peak observed between the HCN-VIB line and the HCO^ 
absorption feature is part of the unabsorbed high velocity HCO^ 
(similar to what is seen in HCN where the absorption is less pro¬ 
nounced). 

In the other galaxies the absorption profiles are instead some¬ 
what blue-shifted by -50 to -100 km s“^ for Zw049.057 and -100 
to -300 km s“^ for Arp220W and IRAS 17208-0014. 

4.2. Vibrationally excited HCN 

We find luminous HCN-VIB J=3-2 or 4-3 If line emission 
in all four sample galaxies. The emission originates from the 
very nucleus of the galaxies with upper limit diameters (FWHM) 
ranging from 37 to 136 pc (Tab. O. For the ULIRG merger 
IRAS 17208-0014 the HCN-VIB emission appears to be asso- 
ciated primari l y with the western of the two nuclei found by 
iMedling et al.l (120141) . For all four galaxies, the HCN-VIB emis¬ 
sion is significantly more compact (by factors of at least 2- 
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Fig. 2. IRAM Plateau de Bure spectra (continuum subtracted) of IC860 (left panel) and Zw049.057 (right panel). The systemic 
velocities of HCN, HCO^ and HCN-VIB J=3-2 are marked with dashed vertical lines. The self- and continuum absorption leads to 
double-peaked HCN and HCO^ y=0 spectra while the HCN-VIB lines are singled peaked close to systemic velocity (see Sect. 14.21) . 
The spectra also show prominent, single peaked lines of HOC^, CH 2 NH, CH 3 OH and HC 3 N. Fluxes are in mJy beam“^ 




Sky Frequency [GHz] Sky Frequency [GHz] 



Sky frequency [GHz] 


Fig. 3. ALMA spectra (continuum subtracted) of Arp220W (left panel) and IRAS 17208-0014 (centre and right panels). The sys¬ 
temic velocities of HCN, HCO^ and HCN-VIB 7=4-3 and the H 2 S lines are marked with dashed vertical lines. The self- and 
continuum absorption leads to double-peaked HCN and HCO^ y=0 spectra while the HCN-VIB and H 2 S lines are singled peaked 
close to systemic velocity (see Sect. 14.21) . Half of the HCN-VIB is blended with the HCO^ line. Note that part of the HCO^ line is 
outside the band for Arp220W. We also detect lines of HC 3 N and HC^^CCN. Fluxes are in mJy beam“^ 


3) than that of HCN and HCO^ and for IC860 the integrated 
strength of the line exceeds that ofHCO^ in the nucleus. 

The HCN-VIB line widths range from 130 to 450 kms“^ 
and the emission peaks close to systemic velocity. This is most 
clearly seen in the two LIRGs where the line blending with 
HCO^ is less significant. The centre (systemic) velocities are 
marked with a dashed line in Figs. |2] and [3] and for IC 860 it 
is estimated to cz=3880+20 kms"^ from the position velocity 
(pV) diagram. The centre velocity for Zw049.057 is very simi¬ 
lar cz=3900+20 km s“^ where the errors also stem from the pV 
diagram (see Sect. 14.41) . For Arp220W the line marks cz=5434 
km s“^ which is 17 km s~^ below th e system ic velocity estimated 
for Arp220W by ISakamoto et al.l (l2008l) . It is suggested that 
the systemic velocity of IRAS 17208-0014 is cz= 12834 kms“^ 
(iDownes et al.l 1 19931) (also marked with a dashed line in Fig. [3]) 
while iGarcia-Burillo et al.l (12015l) suggest that the systemic ve¬ 
locity is blue-shifted by 30 km s“^from this value. We have in¬ 
dicated this lower value of the centre velocity in the pV diagram 
for IRAS 17208-0014 (Fig.©. 

4.3. Other species 

For IC860 and Zw049.057 (see Fig. [2]) we detect the HC 3 N 
7=29-28 y 7=2 line (TEi/k=^^9 K) on the redshifted and the 
(4i 3 - 3 1 , 2 ) transition of CH 2 NH at 266 GHz (with energy level 
TEiik='^2 K) on the blue side of the HCN 3-2 line and sev¬ 


eral blended CH 3 OH lin es between CH 2 NH and HCN-VIB. It 
is interesting to note that lSalter et al.l (l2008l) find bright lines of 
CH 2 NH in their cm-wave spectral scan of Arp220. We also de¬ 
tect very bright HOC^ 3-2 line emission on the blue side of the 
HCO+ 3-2 line. 

For Arp220W and IRAS 17208-0014 (see Fig. [3]) bright 
7=39-38 HC 3 N line emission is blended with the HCN 4-3 line 
and H^^CCCN is found in between the HCN and HCN-VIB line. 
For IRAS 17208-0014 we have also detected luminous line emis¬ 
sion from the 32,1 - 3 1,2 and 43 j - 42,2 transitions of ortho-H 2 S 
and para-H 2 S at 369 GHz. The two lines are shifted only 26 MHz 
from each other (energy levels r£j/j^=137 K and 245 K respec¬ 
tively) so they are effectively blended into one line. 

4.4. Position velocity diagrams 

For IC860 and Zw049.057 the HCN-VIB emission is spatially 
unresolved and peaks at velocities close to systemic (Fig. ID). 
In comparison, the brightest HCN emission occurs outside the 
HCN-VIB at higher absolute velocities, with a velocity differ¬ 
ence between the two peaks of 230-260 kms“^ Even though 
much less pronounced than for the HCN and HCO^ lines, the 
HCN-VIB line also seems to be slightly double-peaked. The ve¬ 
locity shift between the peaks is relatively small: ^^80 km s“^ 

For the ULIRGs Arp220W and IRAS 17208-0014 the HCN- 
VIB line is broad enough for it to be partially blended with 
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Fig. 4. Position velocity (pV) diagrams (continuum subtracted) for IC860 (left) and Zw049.057 (right) where the x-axis shows the 
radial velocity (relative to cz=3880 km s“^and 3900 km s“^respectively) and the y-axis shows the position along the major axis (for 
both galaxies the position angle is -20°). The upper panel shows the HCN, the centre the HCN-VIB and the lower panel the HCO^ 
pV diagrams. The contours are 3.5x(l,3,5,7,9,ll,13) mJy beam“^(IC860) and 3.5x(l,3,5,...17) mJy beam“^ (for Zw049.057). The 
colour ranges from 2 to 30 mJy beam“^ The white vertical dashed line indicate the systemic velocity. The horizontal white line 
indicates the nuclear position of the HCN-VIB line. 
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Fig. 5. Position velocity (pV) diagrams (continuum subtracted) 
for IRAS 17208-0014 where the y-axis shows the position along 
the major axis (position angle 90°) and the x-axis shows 
the radial velocity relative to the cz of 17208-0014 kms“^ 
Upper panel: HCN 4-3. Centre panel: H 2 S Lower panel: 
HCO^ 4-3 blended with HCN-VIB. Colour ranges from 2 to 
70 mJy, contour levels are5x(l,3,5,...17) mJy beam“^ The ver- 
tical dashed line ni a ^rks z= 0.04271 the suggested redshift by 
l&arcia-Burillo et ^ (l2015l) . while zero velocity marks the red- 
shift z=0.04281. The horizontal white line indicates the nuclear 
position of the HCN-VIB line. 


the HCO^ line (see example for IRAS 17208-0014 in Fig. [5]). 
Observations at higher spatial resolution are necessary to com¬ 
pletely disentangle the HCN-VIB emission. The HCN-VIB pV 
diagrams for both Arp220W (to be presented in Martin et al. (in 
prep.)) and IRAS 17208-0014 (Fig. [5]) therefore only show half 
of the line but in both cases a single peak at systemic velocity 
is suggested. For IRAS 17208-0014 the 369 GHz H 2 S line has 
a profile that does not show the pronounced double structure of 
HCN and HCO^. The pV diagram (Fig. [5]) reveals that the H 2 S 
line emission peaks where HCN and HCO^ have a minimum. 
Note also the blue-shifted high-velocity extension to the HCN 
line which may be the HC N counterpart to the CO 2-1 high ve¬ 
locity outflow detected by iGarcia-Burillo et afl (l2015l) . It may 
be similar t o the HCN 3 -2 ou tflow already found in the ULIRG 
Mrk231 by lAalto et alJ (l2015h . 

4.5. Continuum 

Continuum fluxes and FWHM source sizes at 260 or 350 GHz 
can be found in Tab. [3l The FWHM size of the continuum 
emission region is larger than the HCN-VIB emission for two 
sources - Zw049.057 and IRAS 17208-0014. In IC860 the con- 
tinuum source is unr esol ved as is the HCN -VIB emission. 
ISakamoto et all (l2008l) and IWilson et al.l (l20l4 have shown the 
existence of an optically thick and compact submm contin¬ 
uum source in Arp220W consistent with our result. In both 
IRAS 17208-0014 and IC860 the lower limit to the brightness 
temperature of the 260 and 350 GHz continuum sources is 
Tb >40 K suggesting that these sources may also be optically 
thick. 


Table 3. Continuum fluxes and source sizes^. 


Name 

Peak 

[mJy] 

Integrated 

[mJy] 

0 

[pc] 

Tb 

[K] 

IC 860 

35 + 2 

35 + 2 

<37 

>43 

Zw 049.057 

20+1 

44 + 3 

145 X 86 

>5 

Arp220W 

283 + 3 

384 + 4 

62x50 

>160 

IRAS 17208 

87 + 1 

116 + 2 

248 X 162 

>40 


^Continuum levels were determined through a zeroth-order fit 
to line-free channels in the uv-plane. Source sizes (diame¬ 
ters) are full width half maximum (FWHM) two-dimensional 
Gaussian fits. The continuum fluxes and source sizes are de¬ 
termined for 260 GHz for IC860 and Zw049.057 and for 350 
GHz for IRAS 17208-0014 and Arp220. The continuum values 
for Arp220W are fro m Martin et al. (in prep. ) and are consistent 
with those reported in ISakamoto et al.l (120081) . 


5. Discussion 

5.1. Self-absorption in the HCN and HCO^ lines 

The nuclear HCN and HCO^ spectra are double peaked in con¬ 
trast to the HCN-VIB lines (and also the CH 2 NH, H 2 S and 
HOC^ lines). We argue that this is due to a combination of self- 
and continuum absorption in the HCN and HCO^ lines where 
photons at the line centre are reabsorbed by cooler gas in the ex¬ 
ternal parts of the HCN-VIB emitting core and plausibly also in 
more extended regions, (see Fig. [6] for an illustration). A tem¬ 
perature gradient and high column densities will therefore re- 
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Table 2. Line flux densities. 


Line 

Peak^ 

[Jykms“^beam“^] 

Integrated^ 

[Jykms“‘] 

L 

[LoXlO^] 

[pc] 

PA+ 

[”] 

AV 

[kms“‘] 

IC860 

HCN 3-2 

8.8 + 0.3 

19.4 + 0.7 

18.5 

128 X 80 

20 


HCN-VIB 3-2 

3.8 + 0.4 

4.0 + 0.4 

3.8 

<37 


130 

HCO^ 3-2 

3.3 + 0.2 

9.0 + 0.4 

8.6 

133 X 94 

20 


CH 2 NH 

2.2 + 0.2 

3.4 + 0.4 

3.2 

<69 


160 

HOC+ 

1.1 +0.2 

1.3 + 0.3 

1.2 

< 0.2 


140 

Zw049.057 

HCN 3-2 

18.8 + 0.4 

64.4+ 1.5 

60.7 

172 X 137 

25 


HCN-VIB 3-2 

4.2+ 0.2 

6.8 + 0.5 

6.4 

<86 


250 

HCO^ 3-2 

7.3 + 0.2 

25.3 + 0.7 

23.8 

133 X 94 

30 


CH 2 NH 

4.6 + 0.3 

11.3 + 1.1 

10.8 

154 X 80 


250 

HOC+ 

1.7 + 0.2 

2.6 + 0.4 

2.5 

< 0.4 


200 

Arp220W 

HCN 4-3 

130 +5 

294 +15 

670 

310x250 

146 


HCN-VIB 4-3 

55 + 3 

70 +5 

160 

62x50 

133 

350^ 

HCO+ 4-3 

55 + 2 

123 +7 

286 

310x250 

140 


IRAS17208 

HCN 4-3 

53.7 + 0.1 

86.0 + 0.2 

941 

307 X 256 

95 


HCN-VIB 4-3 

8.0 + 0.2 

8.0 + 0.2 

87 

<136 


450^ 

HCO^ 4-3 

30.8 + 0.05 

55.8 + 0.05 

610 

350 X 299 

100 


H 2 S 

36.3 + 0.1 

45.3 + 0.2 

496 

196 X 136 

95 

450 


^Calculated from two-dimensional Gaussian fits to the integrated intensity emission using the AIPS task IMFIT which 
deconvo l ves th e clean beam from the fitted component size. The accuracy and errors of this process are discussed in 

ICondonl i fT^ . 

*From Gaussian fits where the line is blended with HCO^. Limits to HCN-VIB brightness temperatures are LbCHCN- 
VIB)>20 K for IC860, >6 K for Zw049.057, >8 K for IRAS 17208-0014 and 68 K for Arp220W. 

For IC860 and Zw049.057 1"=286 pc, for Arp220W 1"=388 pc and for IRAS 17208-0014 1"=853 pc. 


suit in self-abso rbed lines as seen a lso in embedded sources in 
the Milky Way (iRolffs et al.ll2011bh . The self-absorption oblit¬ 
erates any sign of the hot dense nucleus in the HCN and HCO^ 
lines impeding their use as tracers of the properties of extremely 
dust-embedded n ucle i. The impact of self-absorption is further 
discussed in Sec. l5.2l below. 

How do we know that the twin-peaked HCN and HCO^ lines 
are not caused by a real absence of dense molecular gas in the 
nucleus? We need other lines to compare with to determine the 
cause of the absorption. 

5.1.1. Evidence from the HCN-VIB line 

We can test the self-absorption theory for HCN through compar¬ 
ing the intensity in the vibrational y 2 state (both the le a nd If 
line) to that of the vibrational ground state - 4/7/ in eq. ED. 
Unless the HCN-VIB line is masing, < Ii. The Iu=IQlCN- 
VIB)(l^-l-l/) and for thermal excitation the le line is as bright 
as the If line and hence Iu=2 x /(HCN-VIB). Thus, for lu < Ii 
to hold, the HCN-VIB/HCN line intensity ratio, Kyi\y, should be 
<0.5. Note, that here we are assuming that all of the HCN emis¬ 
sion is emerging from the same gas that produces the HCN-VIB 
line. Any additional HCN-emission would further lower the ’’al¬ 
lowed” !Rvib ratio. 

We can now inspect the nuclear spectra and determine !Rvib on 
a channel-by-channel basis. For IC860 and Arp220W !Rvib ex¬ 


ceeds unity close to systemic velocities and therefore line pho¬ 
tons of both ground state HCN and le HCN-VIB, which appears 
close to the HCN y=0 line centre (see Sect. n, are being ab¬ 
sorbed (see zoomed-in example for IC860 in Fig. [7]). However 
for Zw049.057 the ratio is !Rvib=0-5 and for IRAS17208-0014 
!Rvib=0-25 and additional evidence is required to claim that 
the HCN and HCO^ line shapes are due to self-absorption. 
Furthermore, continuum absorption by the optically thick HCN 
and HCO^ lines can be important, given that the continuum 
brightness is similar to the continuum-subtracted line brightness 
at central velocities. 

5.1.2. Comparing with other lines 

We can use other, optically thin, species to trace the dense gas in 
the nuclear region. In Appendix IA. 1 1 we discuss the CH 2 NH line 
emission in Zw049.057 and the H 2 S line emission of 117208- 
0014 and how they support the notion of self-absorbed HCN 
and HCO^ in these galaxies. In Appendix IA.2I we also discuss 
how the HCO^/HOC^ ratio is also consistent with self-absorbed 
HCO+. 

5.2. Impact of self-absorbed HCN and HCO^ lines in galaxy 
nuclei 

The effect of self absorption on the HCN and HCO^ lines (see 
illustration in Fig. [6]) demonstrates that the emission cannot be 
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viewed as emerging from an ensemble of clouds that are inde¬ 
pendent of each other - but that emission from clouds closer to 
the nucleus is being absorbed by clouds farther out. The notion 
that the spectrum observed towards a galaxy involves ’’counting 
clouds” - that may be individually optically thick but collectively 
’’optically thin” (i.e., not overlapping each other in the 2D (sky 
plane) or the velocity (line-of-sight) space) does not hold for 
these galaxy nuclei. 

This bre aking down of the ens emble approximation has been 
predicted bv iDownes et al.l (1 19931) who suggested that (U)LIRG 
nuclei should be viewed as a single cloud. The ensemble approx¬ 
imation is appropriate when the mean gas surface density is low 
or moderate, but becomes less and less valid when observing a 
dense ULIRG/compact nucleus at increasing resolution. A sim¬ 
ilar conclusion that emission is not emerging from individual, 
no n-overlapping clouds in the nuclei of Arp220 was also drawn 
bv IScoville et al.l (l2015l) . 

When lines of HCN and HCO^ become self-absorbed they can¬ 
not be reliably used to probe the central regions of galaxies on 
scales where the absorption occurs - for our sample galaxies this 
means r < 100 pc. For the studied galaxies this is also where most 
of the emission comes from, in contrast to many nearby galax¬ 
ies with central molecular zones extended over larger regions. 
Thus, the absorbed HCN and HCO^ line profiles are strongly im¬ 
pacting our understanding also of the global properties of these 
galaxies. Very recently self-absor bed CO 6-5 has been d etected 
towards both nuclei of Arp220 bv iRangwala et al.l (l2015l) . 

5.3. HCN-VIB lines as probes of extreme mid-IR surface 
brightness and luminosity 

Vibrationally excited lines of HCN are excellent probes of the 
nuclear dynamics and serve as proxies for the mid-IR luminos¬ 
ity density in deeply obscured nuclei where V(H 2 )> 10 ^"^ cm“^ 
and where dust temperatures exceed 100 K. Mid-IR excitation 
of HCN is favoured in regions of high dust temperature and with 
significant gas and dust column densities. This is because an 
rB(14yum)! 100 K is necessary to start to excite the HCN vibra¬ 
tional ladder and this requires an optically thick mid-IR source 
- which in turn requires an H 2 column density V(H 2 ) ^2 X 10^3 
cm"^ (see Appendix iBl for a discussion). Note that here we as- 
sume that the molecules receive the exciting radiation isotropi¬ 
cally. If they do not then the background source must be even 
hotter to make up for the loss in solid angle. 

The HCN-VIB lines of the sample galaxies are very luminous 
(see discussion in Sect 15.5.11) and lower limits to Tb(HCN-VIB) 
exceed 20 K for IC860 and Arp220W (see footnote to Tab.O. 
Therefore, the filling factor of hot dust and gas is extremely high 
in the inner tens of pc of the sample galaxies - much higher than 
on similar size scales in Milky Way massive star forming regions 
or in extended starbursts. The implication is that the HCN-VIB 
emission we detect towards these galaxies cannot be emerging 
from a collection of Milky-Way-like hot cores. The high filling 
factor of hot, optically thick gas and dust suggests that we can 
use the minimum I4pm brightness temperature of 100 K to esti¬ 
mate a lower limit to the mid-infared luminosity. 

We assume a nuclear spherical ’’cloud” of dust for which 
we adopt a simple blackbody where L=47ir^crT^ (cr is Stefan 
Boltzmann’s constant). Taking r from Tab. [2| (i.e. 6/2) and for 
rB(14yum)=100 K the corresponding L(14pm) is 6 x 10^^ L© for 
IC860, 3 X 10^0 L© for Zw049, 2 x lO^^ L© for Arp220W and 8 x 
10^^ L© for IRAS 17208. In terms of temperature this is a lower 


limit since at 100 K the excitation of the HCN-VIB line will not 
be very efficient. We know from the self-absorption structure of 
HCN and HCO+ (see Sect. IQ) that there is a temperature gradi¬ 
ent and a significant fraction of the HCN-VIB line will emerge 
from regions of higher temperatures even closer to the nucleus. 
Thus, when the radius of the emitting region is known - like it 
is for Arp220W - this value will be a firm lower limit to the 
real luminosity. For Arp220W, high resolution observations the 
submm continuum su ggest an optically thick, peak dust temper- 


ature of 160- 200 K ( 

Downes & Eckart l2007l: 

Sakamoto et al.l 

l2008l;lMatsushita et al. 

20091: IWilson et alJl2014 

) demonstrating 


that an assumption of rB(14yum) of 100 K is too low. We there¬ 
fore conclude that these compact mid-IR nuclei may constitute 
20% - 100% of the total observed L(IR). 


From a lower limit to the rB(14yum) we can estimate an intrinsic 
(unobscured) I4pm surface brightness, E 14 . With the assumption 
of blackbody radiation 

Ei 4 = — = Inr^crT^/nr^ = 2crT^. 


For rB(14yum)=100 K, 1^14=3 x 10^^ L© kpc for a rB(14yum) 
of 150 K 1 : 14 = 1.5 X 10^4 Lq ]^pc-2 fQj. 200 K, 214=5 x 10 ^^ 
L© kpc“^. 

Note again, that a more sophisticated approach to the radiative 
transport and source structure is necessary for a more accurate 
discussion of the surface brightness and mid-IR luminosity. 

5.3.1. Attenuated dust SEDs 

The luminous mid-IR nucleus, which is exciting the HCN-VIB 
line, is hidden inside massive layers of dust and its intense emis¬ 
sion may become absorbed and re-emitted at longer wavelengths 
altering the dust spectral energy distribution (SED). The unatten¬ 
uated SED may be hot (e.g. AGN-like) while the observed SED 
instead appears cooler and more like that of a typical starburst. 

In the conservative (100 K) estimate of L(14yum) above for 
IC860 the unattenuated I4pm continuum fiux from the buried 
source should be ^^2.5 Jy. Howe ver, the observed fi ux density 
is more than a factor 30 lower (iLahuis et al.l l2007h . In Eig. 0 
(lower panel) we show an example of a possible intrinsic SED 
of IC860 together with the observed SED. Another example is 
IRAS 17208-0014 where our HCN-VIB observations imply the 
presence of a ^ 4Jy compact mid-IR source while high resolu¬ 
tion imaging at mid-IR w avelengths only fin ds an extended disk 
of a total fiux of 0.2 Jy (ISoifer et al.ll 2000 h . The authors note 
the striking discrepancy between the 2” mid-IR source and the 
0.”26x0.”32 8.6 GHz radio source which sets IRAS 17208-0014 
aside from the other ULIRGs they studied. The bright HCN-VIB 
line towards the nucleus of IRAS 17208-0014 suggests that the 
solution to this discrepancy is an extreme obscuration of the nu¬ 
cleus at mid-infrared wavelengths. 

The trapped radiation from an embedded source may also 
aid in raising the internal temperature significantly. High column 
densities in all directions lead to the dust being optically thick in 
the infrared and this radiation is unable to escape and diffuses 
outwards by multiple absorption/emission events until the dust 
is optically thin to it s own radiation. This process is discussed in 
iRolffs et M (l2011ah where this ’’greenhouse effect” (Gonzalez- 
Alfonso et al, in prep) is used to explain the high masses of hot 
gas near an embedded protostar. 
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Fig. 6. Upper panel: Illustration of the vibrationally excited 
(HCN-VIB) line emission emerging from a buried mid-IR core 
-surrounded by cooler material giving rise to the self-absorption 
feature in the HCN line. An embedded AGN nucleus would 
be one possible scenario that is powering the bright HCN-VIB 
and causing the temperature gradient with the resulting self- 
absorbed HCN and HCO^ lines. Middle panel: Illustration of 
a possible scenario that may account for the observed HCN 
and HCN-VIB emission in IC860 with a buried compact heat¬ 
ing source, where the assymmetric HCN spectrum is caused 
by an inflo w of gas. The model was generated wi th the tool 
described in iGonzalez-Alfonso & Cernicharol (1 19991). which in - 
cludes line overlaps (iGonzalez-Alfonso & Cernicharol Il997h . 
Note that these examples assume a spherical geometry while our 
results may imply a more torus- or disk-like morphology. Lower 
panel: This is an illustration of how the dust spectral energy dis¬ 
tribution (SED) for IC860 is shifted from shorter to longer wave¬ 
lengths by the absorption of the emission from the hot buried 
source. The dashed curve is the intrinsic unobscured SED and 
the stars and the solid line indicate the observed flux and the 
fitted SED. 


5.4. The nature of the buried source 

The HCN-VIB lines are emerging from highly compact (r < 17- 
75 pc) regions centered on the nucleus of the galaxy. The sphere 



-300 -200 -100 0 100 200 300 

LSR velocity [km s~^] 


Fig. 7. Zoom in on the spectral structure of HOC^, HCO^ (top 
two panels) and HCN-VIB,HCN (lower two panels) of IC860 
showing the self-absorbed line shapes of HCN, HCO^ while the 
other two lines peak at systemic velocity. The intensities have 
been normalized. 


of influence Q of the SMBH in the core of the galaxy is typically 
10-20 pc. Therefore, the HCN-VIB emission is emerging near 
to this region and thus originates in a transition-zone constitut¬ 
ing an possible evolutionary link between the larger scale host 
galaxy and the pc-scale environment near the SMBH. 

A dust-enshrouded AGN would be consistent with the high mid- 
IR surface brightness and a temperature gradient causing the 
HCN and HCO+ self-absorption. The SMBH of IC860 (for 
exaimle) would be accreting at 12 - 30% of its Eddington 
limifl to be the cause of the excitation of the HCN-VIB line 
and 50-100% of the EIR luminosity. An alternative possibil¬ 
ity is that the intense mid-IR emission emerges from a com¬ 
pact starburst. The mid-IR surface brightness of > 5 x 10^^ 
Lo kpc“^ exceeds the limiting Eddington flux of a cool, opti¬ 
cally thick starburst. However, a hot, optically thick starburst 
mode with a bulk temperature T > 200 K may exist accord¬ 
ing to [Andrews & Thompson! (1201 ll) and they propose that its 

Eddington flux is Fem « 

equation (7) where Hg is the gas surface density). 

Eor both cases the HCN-VIB line is tracing a phase of ex¬ 
ceptional growth in the nucleus of the galaxy in the form of a 
rapidly accreting SMBH and/or the fast build-up of a nuclear 
stellar spheriod/disk. 


^ The region where its gravitational potential overwhelms that of the 
galaxy with rg = ICMbh/cts^ (erg is the stellar velocity dispersion). 
It is usually larger than the region inside the Bondi accretion radius 
race = 2 GMbh/cs^ (cs is the speed of sound in the ISM) where the 
gravitational energy of the SMBH dominates over the thermal energy 

of the gas. _ 

^ The K-band luminosity of IC860 dBrown et al.l 120141) implies a 
mass of the SMB H of ^ 2 x 10^ M^ (us ing the prescription described 
in section 3.4 of ISakamoto et al.l (l2013h ). An SMBH accreting at its 
Eddington limit would have L = Lem > M/MqLo. 


9 










































S. Aalto et aL: Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN 


5.4.1. A rotating dust torus? 

Interestingly, the HCN-VIB emission of IC860 and Zw049.057 
also shows a faint signature of a twin velocity peak in the pV 
diagram (see Fig.|4]). This implies that the emission is not arising 
from a sphere or filled disk but is consistent with a torus-like 
structure. However, small (pc-scale) optical depth effects in an 
opaque dust core cannot be excluded, or the possibility that also 
the HCN-VIB line is partially self-absorbed. The width of the 
HCN-VIB line suggests that the enclosed mass inside r=18 pc 
exceeds 5x10^ M© for IC860. The K-band lum inosity implies a 
mass of the SMBH of 2 x 10^ M© (see Sect. 15.4b which has error 
bars of roughly a factor of five. Higher resolution observations 
of the vibrationally excited HCN line will allow determinations 
of the SMBH mass of higher accuracy. 

5.5. Where do we find HCN-VIB emission? 

Even though there has been no systematic survey, the galaxies 
with HCN-VIB emission detections so far are either ULIRGs or 
LIRG early type spirals (see Tab|4]). 

The detection of HCN-VIB emission requires the presence 
of high masses of gas and dust and a compact heating source. 
Early type spirals have steeper potential wells than late type spi¬ 
rals which may more easily facilitate the formation of compact 
dust cores. However, optically thick, dusty nuclei should be eas¬ 
ier to find in early type spirals than in late type ones since the 
former generally do not have contaminating star formation in 
their disks. A survey will reveal whether the HCN-VIB lines 
preference for early type spirals is real o r a result of observa- 
tional bias, but it is interesting to note that iMignoli et al.l (l2Q13l) 
find that moderate luminosity, obscured type-2 AGNs occur pri¬ 
marily in early-type Sa-Sb galaxies. 

The link between major mergers and ULIRGs is well known 
( Sanders & Mirabellll99d iBekki et al.lll999l: lElbaz & Cesarskvl 

l2003h and luminous, dense circumnuclear disks of gas and stars 
are a natural consequence of g as-rich galaxy mergers since the y 
trigger massive inflows of gas (lBarnestt2002l:lMaver et al.ll2007l) . 
In the very inner region the gas may contribute to the rapid 
growth of SMBHs and/or the build up of nuclear stellar popula¬ 
tions. The incidence of HCN-VIB line emission in ULIRG ma¬ 
jor mergers will give important clues to when during the ULIRG 
phase the nuclear region is undergoing accelerated growth. 

5.5.1. The L(HCN-VIB)/L(FIR) ratio 

In Tab. |4] we list the L(HCN-VIB)/L(EIR) ratio for the eight 
galaxies that have HCN-VIB detections so far: four lenticu- 
lars (NGC1377, NGC4418, IC860 and Zw049.057) and four 
ULIRGs (IRAS20551-4250, Arp220, Mrk231, IRAS17208- 
0014). We have also added the Galactic massive star formation 
region SgrB2(M) to illustrate that the ratios for most of the HCN- 
VIB detected galaxies are higher by an order of magnitude than 
those of a forming stellar cluster in the centre of the Milky Way. 
Even if we populated the inner region of (for example) IC860 
solely with SgrB2(M)-like regions this would not reproduce our 
observations. We have also added the ULIRG IRAS 13120-5453 
to Table |4] since we have high signal-to-noise ALMA observa¬ 
tions (Privon et al in prep.) with a good upper limit to the HCN- 
VIB line. 

There is no obvious relation between L(HCN-VIB)/L(LIR) and 
LIR luminosity. The ratio will naturally be dependent on the tem¬ 
perature of the exciting dust source and the fraction of the L(LIR) 


emission that is emerging from the same region as the HCN-VIB 
emission. However, a trend to be explored further is the cor¬ 
relation with the strength of molecular outfiows (see Table ID). 
The rapidly evolving nature of the studied nuclei is illustrated 
by what seems to be in- and outfiow signatures in the self- 
absorbed HCN and HCO^ lines, but also in other absorption- 
and emission lines. All galaxies with a detected HCN-VIB line 
show evidence of non-circular motions: inflows, outflows or 
both. Tentatively, the galaxies with the relatively brightest HCN- 
VIB lines show signs of slow non-circular motions: NGC4418, 
Arp220W, IC860, Zw049.057, IRAS 17208-0014 (see Tab.^for 
references). (The reversed P-Cygni profiles we find for IC860 
could be evidence of infall, but more observations and modelling 
are required to determine their precise nature.) Galaxies with fast 
outfiows (i.e. where the outfiow velocities exceed the galaxy’s 
rotational (virial) velocity) have fainter HCN-VIB signatures. 
The ULIRGs Mrk231, IRAS20551 and IRAS13120 have promi¬ 
nent molecular outfiows, seen in the OH 119 jam line, with max¬ 
imum velocities >1000 kms“^(see Tab. |4] for references). The 
lenticular galaxy NGC1377 has a prominent molecular outfiow 
in relation to its relatively modest escape velocity and mass 
(lAalto et alJl2ni 2h . 

Apparently in contrast, ISooon et^ (l2013h find, using the 
119 yum OH absorption line, that the highest outfiow velocities 
are found in galaxies with deep mid-IR silicate absorption indi¬ 
cating that they are embedded. It is possible that the HCN-VIB 
line is probing a deeper, more extreme form of nuclear obscura¬ 
tion - and possibly also more transient - than the one trace d by 
the mid-IR silicate absorption. Indeed, iRoche et ^ (1201 5l) find 
for NGC4418 that a large fraction of the silicate absorption is 
arising f rom ’’cool material outside the w arm molecular core”. In 
addition JConzalez-Alfonso et al.l (1201 5l) find that the silicate ab¬ 
sorption is biased towards relatively unabsorbed mid-IR emitting 
regions and propose high-lying OH absorption lines as tracers of 
embeeded, hot nuclei. (The correlation between HCN-VIB and 
these highly excited OH lines will be explored in future work.) 

A potential scenario is therefore that the HCN-VIB phase 
peaks right before the nuclear feedback mechanism is fully 
developed. This would be consistent with our conclusion in 
Sect. 15.41 that the HCN-VIB line is probing a transient phase of 
rapid nuclear growth. A possible complication to this picture is 
the recent discovery of a 800 km s“^ CO 2-1 outfiow in the HCN- 
yiB luminous ULIRG IRAS 17208-0014 dGarcia-Burillo et all 
l2015l) . No fast 119 fim OH wind has been found and it is possi- 
ble that these absorption studies may miss in- and outfiows that 
occur deep inside an optically thick dust shell. 

5.5.2. HCN-VIB as a survey tool for obscured high-surface 
brightness activity 

Our observations show that luminous extragalactic HCN-VIB 
emission is emerging from compact nuclei of high mid-IR sur¬ 
face brightness and that vibrationally excited HCN is a strong 
spectral signature of a nucleus in rapid growth - either through 
Eddington limited (or near-Eddington) nuclear starburst activity 
or of accreting AGNs. The HCN-VIB line is an effective survey 
tool to find high surface brightness nuclei hidden behind high 
masses of cooler dust and gas - even with observations at low 
spatial resolution. Detecting the vibrationally excited HCN also 
allows us to reconstruct the intrinsic, hotter SED of the luminos¬ 
ity source. 

The HCN-VIB line from a nuclear compact core will show 
up even if there is more extended star formation which masks 
and confuses other tracers such as the [C II]/LIR ratio and radio 
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Table 4. L(HCN-VIB) vs L(EIR) 



Name 

log L(FIR) 

L(HCN-VIB) 

L(FIR) 

[10-8] 

Molecular 

Outflow^ 

Ref 

7=3-2 

NGC4418 

11.08 

3.8 


[1] 

IC 860 

11.17 

3.2 


This work 

Zw 049.057 

11.27 

3.4 


This work 

Mrk231 

12.37 

0.5 

Y 

[2] 

SgrB2(M) 

6.8 

0.3 

y 

[3] 

7= 4-3 

NGC1377 

10.13 

1.0 

Y 

[5] 

NGC4418 

11.08 

9.5 

n* 

[1] 

IRAS20551 

12.00 

0.4 

Y 

[4] 

Arp220W 

12.00 

15 

y 

This work 

IRAS 13120 

12.26 

<0.2 

Y 

[6] 

IRAS 17208 

12.39 

3.6 

y 

This work 


rn ISakamoto et alJ llM(l l2013h. [21 lAalto et alJ ( |2015h . [3] 
iRolffs et al.N20TM . [41 llmanishi & Nakanishil <2^ . [5] Aalto et 
al (in prep.), [6] Privon et al. (in prep.) 

capital Y indicates molecular outflow velocities exceeding the 
escape velocity for the nuclear re gion of the galaxy. References 
for the outflow s/ infall: NGC4418 (iGonzalez-Alfonso et al.l 120121: 
Sakamoto et al.1 120131 : ICostagliola et al.1 [20131: IVeilleux et^ 
2013h : IC869 - this work; Zw049 .057 iFalstad et al.l 1201^ 
IRAS 17208-0014 ([yeillei^ et al.ll2013h (but note recent work by 
iGarcia-Burillo et al.l (120151) wh ere CO 2-1 outflow velocities of 
800 I on s~^are found), Mrk231 jFischer et al.ll20TQl : IVeilleux et al.l 
120131 : IGonzalez-Alfonso et al.l |2014|) , IRAS 13120-5453 and 
IRAS20551 (IVeilleux et al.1 l2013h . ArD220W j^kamoto et al.l 
120091 : IVeilleux et al.ll2013h . SgrB2(m) jRolffs et al.ll2011hh ~ 

^Minor axis dust lanes in the form of cones suggesting current or 
past outflows have been found in NGC4418 (ISakamoto et al.|[20l^ 
and also in IC860 (unpublished Nordic Optical Telescope U,B 
images by F. Costagliola). Also in Zw04 9.057 a peculiar mino r axis 
outflow-like dust feature has been found lScoville et all (l2000h . 


continuum imaging. Galaxies with HCN-VIB emission there¬ 
fore constitute a sample where to search for deeply obscured 
AGN activity. Any X-ray emission should be severely attenu¬ 
ated - even hard(> 5 keV) X-ray emission if column densities 
exceed A^(H 2 )> 10^"^ cm“^. Follow-up, high resolution HCN- 
VIB and continuum studies will help measure the luminosity 
density and enclosed mass which will put strong constraints on 
the nature of the buried source. Thus, studies of vibrationally 
excited lines can be used to reveal and probe a population of 
deeply dust-buried AGNs promising to shed light on how the 
SMBH growth connects to that of the stellar bulge and to enable 
studies of the competition between star formation and accretion. 
Emerging new techniques utilizing He^ recombination lines may 
further ass ist in separating accreting SM BHs from hot compact 
starbursts (IScoville & Murchikovall2013l) . 

6. Conclusions 

- High resolution (0.”4) IRAM PdBI and ALMA mm and 
submm observations of the LIRGs IC860 and Zw049.057, 
and the ULIRGs Arp220 and IRAS 17208-0014 reveal lu¬ 
minous vibrationally excited HCN J=3-2 and 4-3 V 2 =lf 
(HCN-VIB) line emission emerging from buried, compact 
(r < 17-70 pc) nuclei. They have very high inferred mid- 


infrared surface brightness >5x10^ Lopc“^ and are powered 
by dust enshrouded accreting SMBHs and/or hot (> 200 K) 
extreme starbursts. The HCN-VIB emission may be the sign¬ 
post of an enshrouded nucleus undergoing rapid evolution 
either through the growth of its SMBH or a compact stellar 
spheroid. In either case, the HCN-VIB emission is emerg¬ 
ing from inside the inner few tens of pc which consitutes an 
important transition region and a potential evolutionary link 
between the larger scale host galaxy disk and the pc-scale 
environment near a SMBH. 

- In contrast, we show evidence that the ground vibrational 
state rotational lines of HCN and HCO^ J=3-2 and 4-3 
(standard tracers of dense gas in galaxies) fail to probe the 
highly enshrouded, compact nuclear regions due to strong 
effects of self-absorption (caused by a temperature gradient 
and large line-of-sight column densities) and absorption of 
continuum. We advocate HCN-VIB lines as better tracers 
of the dynamics, mass and physical conditions of (U)LIRG 
nuclei when H 2 column densities exceed 10 ^^ cm“^ and 
when there is a buried compact luminosity source. The vi¬ 
brationally excited HCN allows us to identify galaxies where 
nuclear mid-infrared emission may be strongly absorbed and 
re-emitted at longer wavelengths - and to reconstruct the in¬ 
trinsic, hotter dust spectral energy distribution (SED) of the 
buried luminosity source. 

- With this study HCN-VIB lines have now been found in 
eight (U)LIRGs with AGN signatures and/or evidence of 
compact mm/submm continuum emission. They are either 
ULIRG mergers or early-type spirals with centrally concen¬ 
trated molecular gas and dust. The embedded galaxy nuclei 
show strong signatures of non circular motions that could 
be interpreted as infall and/or molecular outflows - also in 
the self-absorbed HCN and HCO^ spectra. IC860 have re¬ 
versed P-Cygni profiles that may be evidence of inflowing 
gas. A tentative conclusion is that galaxies without evidence 
of fast outflows (i.e. velocities greater than rotational) have 
the strongest HCN-VIB line in relation to their IR luminosi¬ 
ties. This implies that the HCN-VIB emission peaks before 
the onset of strong feedback. 

- All galaxies with detected HCN-VIB lines show very rich 
molecular spectra and we report the detection of luminous 
emission lines of CH 2 NH and HOC^ and vibrationally ex¬ 
cited HC 3 N in IC860 and Zw049.057 also lines of CH 3 OH 
are detected. In IRAS 17208-0014, very bright H 2 S emission 
is found. 
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Appendix A: Self-absorbed HCN and HCO^ lines: 
Comparing with other species 

A.1. Comparing with the H 2 S and CH 2 NH lines 

Both the CH 2 NH and H 2 S molecules have high dipole mo¬ 
ments (/ia=1.340 D and /j=0.91 D respectively) and the transi- 
tions observed here have critical densities i n excess of 10^ cm“^ 
(lDickensetalJll997l: ICrockett et al']l20l4 The CH 7 NR/HCN 
abundance ratio is typically 10"^ - 10"^ dPickens et alJll99% 


suggesting that the CH 2 NH line emission should suffer self¬ 
absorption only in extreme conditions. We find that the CH 2 NH 
line in Zw049.057 has a single peak at systemic velocity and 
the emission has a maximum at the same spatial position as the 
HCN-VIB line. Similarly, in IRAS 17208-0014 the 369 GHz H 2 S 
line clearly peaks at systemic velocities and /(H 2 S)>/(HCN) 
between 0 and -200 kms“^ Even though it is possible for the 
abundances of the t wo molecules to be si milar in the inner re¬ 
gions of hot cores (ICrockett et al.l l2014l) the H 2 S abundance 
must be two orders of magnitude greater than that of HCN 
for /(H 2 S)>/(HCN). This would constitute an hitherto unknown 
astrochemical scenario. Thus the simplest explanation is that 
HCN is self-absorbed in the inner region of Zw049.057 and 
IRAS 17208-0014. 

A.2. Comparing with the HOC^ line 

Zw049.057 and IC860 show intensely bright HOC^ 3-2 line 
emission with /(HOC^)>/(HCO^) at systemic velocity (see ex¬ 
ample for IC860 in Fig. [7]). Assuming the excitation of the two 
molecules is similar, this either requires that HOC^ abundances 
exceed those of HCO^ - or that HCO^ is self absorbed. The 
abundance ratio between the two molecules is a measure of the 
electron abundance and in Galactic dense clouds HCO^/HOC^ 
is typically a few times 10^. In more extreme regions such as 
the photon dominated regions (PDRs) of the starburst galaxy 
M82 and near the nucleus of the Seyfert gal axy NGC1068, this 
ratio is found to range between 44 and 1 (iFuente et al.l l2008l: 
lUsero et ^l2004l) . In no source has HOC^ yet been found to be 
more abundant than HCO^ and chemical m odels do not find this 
even i n extreme X-ray d ominated regions (ISpaans & Meijerink 
I 2 OO 7 I: iBavet et ^|2011|) . We therefore conclude that HCO^ is 
self-absorbed in Zw049.057 and in IC860. 

Appendix B: Why does HCN-VIB thrive in obscured 
regions? 

The HCN y 2 =l vibratio nal ladder can commence being excited 
when rB(14/im)! 100 K (ICarroll & Goldsmithlll98ll:lAalto et al.l 
l2007h . The mid-IR emission must therefore be opaque when the 
dust temperature is Tdust^lOO K. This requires an optical depth 
of at least T(14yum)^5 and for a standard Galactic extinction 
curve (IWeingartner & Drainell200ll) and Ay ~ Mj/2x 10^^ cm“^ 
this is equivalent to an H 2 column density A^(H 2 ) ^2x10^^ cm"^. 
For higher Tdust a lower r will be sufficient for rB(IR)=100 K to 
be satisfied, but at a dust temperature of 400 K an H 2 column 
density in excess of = 10^^ cm“^ is still necessary. 

The vibrational temperature Tvib is a measure of the HCN- 
VIB excitation temperature and hence of the brightness of the 
background mid-IR radiation field the molecules ’’see”. We can 
estimate an effective Tvib through comparing the intensities of 
the HCN-VIB (4) and vibrational ground state lines (//) assum¬ 
ing optically thin emission, 

4/7/ = AulAi(Vulvifgulgi&w(-EulTvib)- (B.l) 

For a Tvib of 100 K the optically thin intensity ratio between 
the HCN-VIB and the HCN line is 3 x 10”^ and the HCN-VIB 
line is not detectable. Thus, high column densities are required 
to observe the HCN-VIB line when Tyib^lOO K. In contrast, for 
Tvib=400 K the optically thin ratio between the HCN and the 
HCN-VIB line is 0.08. However, such a high Tb (mid-IR) re¬ 
quires rdust=400 K combined with a large dust column density 
for optically thick emission - or a very large Tdust > 400^. Hence 
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mid-IR excitation of HCN is favoured in regions of high dust 
temperature and with significant gas and dust column densities. 
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